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ABSTRACT 

Aims. Angular momentum (J) loss requires magnetic interaction between the forming star and both the circumstellar disk and the 
magnetically driven outflows. In order to test these predictions many authors have investigated a rotation-disk connection in pre-main 
sequence objects with masses larger than about 0.4 M Q . For brown dwarfs (BDs) this connection was not investigated as yet because 
there are very few samples available. We aim to extend this investigation well down into the substellar regime for our large sample of 
ss 80 BDs in the Orion Nebula Cluster, for which we have recently measured rotational periods. 

Methods. In order to investigate a rotation-disk correlation, we derived near-infrared (NIR) excesses for a sample of 732 periodic 
variables in the Orion Nebula Cluster with masses ranging between w 1.5-0.02 M Q and whose IJHK colors are available. Circumstellar 
NIR excesses were derived from the A[I-K] index. We performed our analysis in three mass bins (>0.4M G , 0.4-0.075 M B , and 
< 0.075 M G ). 

Results. We found a rotation-disk correlation in the high and intermediate mass regime, in which objects with NIR excess tend to 
rotate slower than objects without NIR excess. Interestingly, we found no correlation in the substellar regime. A tight correlation 
between the peak-to-peak (ptp) amplitude of the rotational modulation and the NIR excess was found however for all objects with 
available ptp values (< 0.4 M Q ). We discuss possible scenarios which may explain the lack of rotation-disk connection in the substellar 
mass regime. One possible reason could be the strong dependence of the mass accretion rate M on stellar mass in the investigated 
mass range (M a M 2 m 28 ), which is expected to result in a corresponding mass dependence of J/ J. 

Key words, stars: low-mass, brown dwarfs, pre-main sequence - stars: rotation, starspots - technique: photometric 



1. Introduction 

At present observations and theory together give a comprehen- 
sive picture of star formation. Nevertheless, how young pre-main 
sequence (PMS) objects effectively lose angular momentum is 
still not fully understood. PMS stars increase their mass by ac- 
creting material from the circumstellar disk. At this stage, accre- 
tion and ongoing contraction would contribute to a net spinning 
up of the PMS star, which is expecting to rotate at nearly breakup 
velocity (e.g. Vogel & Kuhi, 1981). Yet, numerous observations 
reveal that these objects rotate much slower, with typical rota- 
tional velocities of ^10% the breakup velocity (e.g. Herbst et 
, al. 120071 and references therein). It is clear that considerable 
amounts of angular momentum are somehow removed from the 
star/disk system during the accretion phase. 
Various theories, all of them requiring ongoing mass accre- 
tion from the disk to the star in combination with magnetically 
driven outflows, are usually proposed to explain this enormous 
removal of angular momentum and the large spread of rotational 
velocities observed (e.g. Hillenbrand 1997, Herbst et al. 2002 
hereafter H2 002; Lamm et al. 120041 EES! Scholz & Eisloffel 
120041120051 Rodriguez-Ledesma et al. |2009a| hereafter Paper 1). 
The proposed theories rely on a disk-locking mechanism (e.g. 



* Full Table 2 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



Camenzind fT990l K onigl fTWTT ). X-winds (e.g. Shu et al. [19941 
Ferreira et al. 2006 ), disk winds or massive stellar winds work- 
ing in conjunction with magnetospheric accretion (e.g. Matt 
& Pudritz, I2005bl ) to account for the angular momentum loss. 
Irrespective of the specific model, there is a general consensus 
that circumstellar accretion disks are somehow responsible for 
the removal of angular momentum, and therefore accreting stars 
should on average rotate slower than non-accreting ones. 
There is wide observational evidence that most low mass pre- 
main sequence stars and even substellar objects have circum- 
stellar disks (e .g. Hi llenbrand et al. [19981 Muench et al. 120071 
Mohanty et al. 120041 Rebull et al. 120061 etc). Accretion features 
have also been observed in substellar objects (e.g. Mohanty et 
al. 120041 1. 

For more than 15 years, several authors investigated possible 
correlations between rotation rates and the presence of circum- 
stellar accretion disks (e.g. Edwards et al. 1993; Rebull, 2001; 
H200 2; Lamm et al. 120031 Rebull et al. 120061 Cieza and Baliber, 
120071 ). While many of them found a rotation-disk correlation for 
solar and low mass stars (e.g. H2002; Lamm et al. (2005 ; Rebull 
et al. 120061 Cie za and Baliber, [2 0071 B iazzo, 120091 . others did 
not (e.g. Rebull, [2001] Ma kidon [20041 Nguyen et al. 120091 . As 
stated by Cieza & Baliber (2007 ), intrinsic dependence of rota- 
tion on mass and age, which were not taken into account in many 
of the previous works, were probably one of the reasons for the 
discrepant results. An additional reason is probably the use of 
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different and sometimes ambiguous accretion indicators. Most 
of the above mentioned studies are using near-infrared (NIR) ex- 
cess as a circumstellar disk indicator. While NIR excess probes 
the presence of circumstellar dust, it does not probe ongoing ac- 
cretion. However, a remarkable correlation between mass accre- 
tion indicators such as Ha emission line widths and NIR ex- 
cess has been reported in the literature (e.g. Sicilia-Aguilar et al. 
120051 l2006l l and therefore using NIR excess as a mass accretion 
indicator is a reasonable approach. 

Muench et al. ( 2001 ) found that about 50% of their brown dwarf 
candidates in the Trapezium region show near-infrared (NIR) 
excesses indicative of a circumstellar disk. Lada (2000) found 
that ^50% of objects in the Trapezium cluster have a JHK ex- 
cess, while their observations at ^3/zm revealed a much higher 
disk fraction of about 85%. Since longer wavelengths provide a 
better disk indication than observations at NIR wavelengths, the 
use of Spitzer data (IRAC bands: 3.6/zm, 4.5 //m, 5.8 /mi, and 
8 fim) seems to be the solution towards an unambiguous disk in- 
dicator. However, the evaluation of the IRAC data of the ONC 
encounters two serious problems, namely a very strong nebu- 
lar background and high stellar densities. These two problems 
obviously get stronger with decreasing distance to the cluster 
center. As stated by Rebull et al. (2006]!, due to the strong neb- 
ular background and source confussion, even relatively bright 
periodic variables close to the Trapezium are missing IRAC de- 
tections. Megeath et al. (2005) estimated that about 700 stars 
may be missing IRAC detection in the inner 5x5 arcmin 2 center 
of the ONC. Only objects in the ONC with spectral types ear- 
lier than M2 (i.e. M ^ 0.4 M if Baraffe et al. (119981 1 models are 
used) were analyzed and searched for the presence of a disk and 
a rotational-disk connection by means of IRAC data (Rebull et 
al. 12006 and Cieza & Baliber [20071 ). Since our sample of peri- 
odic variables is three magnitudes deeper than previous studies 
and contains many objects in the inner most part of the ONC, we 
have to rely on NIR data to identify the circumstellar disks in our 
sample. Although NIR observations are less efficient in detect- 
ing disks than observations at longer wavelengths, the relatively 
high disk fraction of »50% found by Lada (2000 ) and Muench 
et al. (1200 11 1 among large samples of low mass stars and brown 
dwarfs by means of JHK excesses provides valid support that 
NIR excess is a suitable disk indicator for low mass and substel- 
lar objects. 

A considerable number of studies, which have explored the con- 
nection between the presence of disks and rotation in low and 
very low mass stars, are summarized in Herbst et al. 2007] and 
Bouvier, 120071 The major aim of this study is to constrain for 
the first time a possible rotation-disk connection in the substel- 
lar mass regime, and compare these results with higher mass ob- 
jects. This study would not have been possible without the rota- 
tional period measurements of 124 BD candidates in Paper 1. 
In Paper 1 we confirmed the strong period dependence on stel- 
lar mass found for higher mass stars by H2002 and extended it 
into the brown dwarf regime down to as 20M JU/; . We found that 
substellar objects tend to rotate faster than their stellar counter- 
parts. In order to constrain a possible rotation-disk connection, 
in particular for the so far very poorly investigated very low mass 
and substellar regime, we analyze here a sample of 732 periodic 
variables from Paper 1 and Herbst et al. 120021 which have avail- 
able JHK colors to be used as disk indicators. In Sect. 2 we give 
details of the used NIR data sets. Section 3 describes how we 
estimated the masses in our sample. Section 4 defines the pro- 
cedure to account for objects with and without NIR excess. We 
present our results on the correlations between rotational peri- 
ods and peak-to-peak amplitudes with NIR excess in Sect. 5 and 
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Fig. 1. Spatial distribution of all 732 periodic variables with 
near-infrared counterpart. Black squares are 113 objects with 
VLT JHK measurements, red crosses are the 210 periodic vari- 
ables with UKIDSS JHK measurements and blue open circles 
are the 409 periodic variables with 2MASS measurements. The 
circle centered on 1 Ori with a radius of ~ 6!7 represents the 
cluster radius (see Paper 1 for details). 

6 and briefly address the disk fraction found in Sect. 7. In Sect. 
8, we discuss in detail possible scenarios which can result in a 
lack of correlation between NIR and rotation and summarize our 
results in Sect. 9. 

2. Optical and near-infrared data 

The I band data were obtained with the wide field imager (WFI) 
at the MPG/ESO 2.2-m telescope on La Silla, Chile (see Paper 1 
for details). We used the combined (i.e. Paper 1 and H2002) I 
band data, which include 746 periodic variables with masses 
ranging from about 1 .5 M (H2002) down into the BD regime. 
We combined these optical data with JHK data from three dif- 
ferent sources: UKIDSSfl VLT (MacCaughrean et al. in prep.), 
and 2MASS. Ideally we would have used only UKIDSS since it 
provides much deeper and much more accurate JHK data than 
2MASS, but unfortunately only less than 50% of our ONC WFI 
field (33 x 34 arcmin 2 , center approximately at the position of 1 
Ori) is covered in all three bands by the UKIDSS survey. In par- 
ticular the innermost part of the ONC is not covered (for details 
see Fig. 1 ). We therefore used for a fraction of the inner ONC re- 
gion deep VLT-ISAAC data (MacCaughrean et al. in prep.) cov- 
ering approximately the central 7x7 arcmin 2 . For those regions 
in our WFI field which were neither covered by UKIDSS nor by 
the VLT data we used 2MASS data. The cross-identification be- 
tween the periodic variables and these near-infrared catalogues 
was done by looking for positional matches within a 2"search 
radius. Finally 732 out of the 746 known periodic variables from 
Paper 1 have an infrared counterpart for which data in J, H, and 
K are available. 

1 The UKIDSS project as defined in Lawrence et al 2007 uses the 
UKIRT Wide Field Camera (WFCAM; Casali et al. 120071 . 
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Fig. 2. Histograms showing the magnitude distribution of the periodic variables in I, J, H, and K. The three different NIR data sets 
finally used in this paper are shown separately. The I band data from Paper 1 are shown for comparison purposes. 



Many objects have JHK data in more than one of the used cata- 
logues. Due to their higher accuracy we gave higher priority to 
the UKIDSS and VLT data than to the 2MASS data (see Sect. 
2.2 for details). 

The basis of the VLT catalogue are VLT-ISAAC images. This 
catalogue is supplemented and calibrated with data from a size- 
able number of other surveys which have been transformed to the 
2MASS system resulting in a homogeneous catalogue. To allow 
a proper comparison, we transformed the UKIDSS magnitudes 
into the 2MASS system using the relations described in Hewett 
et al. (120061 1. The magnitude distribution of the periodic vari- 
ables in I and JHK are shown in Fig. 2 with separate distributions 
for the three infrared data sets used in this paper. As summarized 
in Table 1, we finally used 210, 113, and 409 JHK measurements 
from UKIDSS, VLT and 2MASS, respectively. Most objects in 
the innermost part of the ONC, for which 2MASS data had to 
be used, are bright objects from H2002, 95% of which have J 
magnitudes brighter than 14 mag. 

2. 1 . Photometric errors of I-K, J-H and H-K colors 

As we will outline in Sect. 4, our IR excess determinations are 
based on I-K, J-H, and H-K colors. Therefore we have to briefly 
discuss the accuracy of these color determinations. 
There are basically two sources of errors: photometric errors of 
the individual I, J, H, and K measurements and errors resulting 



from variability of the periodic variables in combination with the 
non-simultaneous data acquisition in the used bands. 
Non-simultaneous data acquisition occurred for all I-K colors. 
As outlined in Paper 1 (see there Fig. 14), the peak-to-peak (ptp) 
amplitudes in the I band for periodic variables of 0.075-0.4 M Q 
have a median value of about 0.06 mag, and for 85% of them the 
ptp values were below 0.2 mag. Since we will rarely observe the 
periodic variables in minimum in one band and in maximum in 
the other band, I-K colors are expected to have median errors due 
to variability of at most 0.06 mag. For the J-H and H-K colors, 
it is irrelevant in the 2MASS data, which are practically taken 
simultaneously. This means that only for the UKIDSS and VLT- 
ISAAC data one has to investigate the relevance of this effect. 
Since the ptp amplitudes are expected to be significantly lower 
in the NIR than in the I band, the J-H and H-K colors calculated 
from the UKIDSS and VLT data are expected to have median 
errors due to variability of < 0.05 mag. 

The expected mean photometric errors depend on the sur- 
vey/catalogue used. I band photometry from Paper 1 has typical 
photometric errors for a 1=14, 1=18 and 1=20 of m 0.001 mag, 
ss 0.01 mag, and « 0.03 mag, respectively. 
The UKIDSS data are very deep and accurate with typical pho- 
tometric errors for a J=14 , H=13.5 and/or K=13mag object of 
~ 0.002 mag and typical errors of =s0.01 mag for 2 mag fainter 
objects. As stated above, observations in the three bands are not 
taken simultaneously. Nevertheless, as already discussed, we do 
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Table 1. Number of periodic variables from Paper 1 with a near- 
infrared counterpart 





Periodic variables 




At least one band 


All JHK bands 


UKIDSS 


419 


210 


VLT 


113 


113 


2MASS 


684 


409^, 



" 684 periodic variables in our sample have 2MASS JHK colors, but 
since UKIDSS and VLT data have higher priority than 2MASS, 2MASS 
colors have been used only for the 409 objects without UKIDSS and 
VLT data. 

not expect variability to cause errors in extent of 0.05 mag. 
The 2MASS data have photometric errors which are on aver- 
age one order of magnitude higher than those from UKIDSS for 
the same magnitude. As already stated above, objects for which 
2MASS data had to be used are those located in regions of our 
WFI field that are not covered by neither UKIDSS nor VLT (see 
Fig. 1). 

The VLT catalogue reaches 5cr point source limiting magnitudes 
at approximately 22, 21, and 20 for J,H, and K, respectively. 
This is about 2 mag fainter than UKIDSS, and typical photomet- 
ric errors are therefore expected to be «5 times smaller than 
those of UKIDSS for the same magnitude. As already men- 
tioned the VLT data were not taken simultaneously. However, 
each data point is the average of several observations separated 
over timescales of hours to years, and consequently "variability 
errors" will be smoothed out to some degree. 
We expect uncertainties in the I-K color determination domi- 
nated by the non-simultaneous data acquisition of < 10%. From 
the typical photometric errors we expect uncertainties in the J-H 
and H-K colors calculation of about 1 % for the VLT data and 
for the 2MASS data between 3% and 10% for a K= 12 mag 
and K= 14 mag objects respectively. The colors derived from 
UKIDSS data are affected by both photometric and variability 
errors. The J-H and H-K color errors are expected to be clearly 
smaller than 5% and 7% for a K = 13 and 15 mag object, re- 
spectively. This means that in principle we are able to detect 
excesses from the intrinsic J-H and H-K colors which are larger 
than about 1-5% for objects with masses of « 0.1 M and 7-10% 
for the less massive objects in our sample. A rough estimate of 
the amount of NIR excess under typical stellar + disk system pa- 
rameters (e.g. Meyer et al. |1997l and Hillenbrand et al. 1 19981 1 re- 
sults in H and K band excesses which are expected to be of about 
20%. Although the amount of NIR excess is highly dependent on 
the different stellar + disk parameters, the estimated errors in the 
colors should not significantly affect our disk detectability. 

3. Mass estimates: Color-magnitude diagrams 

To estimate masses from evolutionary models we need to know 
individual extinction values, which requires the knowledge of 
the spectral type. Since the spectral types of most of our periodic 
variables are unknown, individual extinction values cannot be 
determined. For the sake of simplicity we assumed in Paper 1 a 
uniform extinction value of A,. = 1.4 mag for all stars in the field, 
and we estimated masses based on the I band values. Although 
we still lack information on spectral types and therefore indi- 
vidual extinction values of our sample, we can make use of the 
near-infrared colors to perform additional/independent mass es- 
timates, which are probably more accurate than the estimates 



based on the I band alone, due to the smaller extinction in the 
near-infrared bands and the smaller variability in this wavelength 
range. 

We used the J vs I-J color-magnitude diagram for the mass esti- 
mates. In principle, H and K bands would be even less affected 
by extinction than the J band, but in these bands infrared excess 
due to stellar disks can pretend a much higher stellar flux than 
is actually present. We investigated the median extinction val- 
ues for the inside and outside R c iuster in the literature (see Fig.Q] 
and Paper 1). First, we computed the median values of the ex- 
tinction provided by Hillenbrand ( 1997 ) in both spatial regions, 
which resulted in A v =1.3 mag and A,,=0.7 mag for the inner and 
outer regions, respectively. Second, we checked extinction maps 
by Chaisson & Dopita (119771 1. The latter reveals that there is no 
strong increase in the extinction towards the inner regions, with 
median extinction values of about 1 and 1.5 in the innermost 
parts of the ONC and below 1 in the outer regions. On the basis 
of these studies we assumed an average A v of 1 .4 mag for objects 
located inside R c iuster and 0.7 for objects outside Rciuster- We used 
the assumed extinction values together with a 1 Myr isochrone 
by Baraffe et al. (19980 to select objects with masses higher 
than 0.4 M (i.e. »M2 or earlier), between 0.075-0.4 M Q (i.e. » 
M2 to M6), and < 0.075 M Q (i.e. later than * M6). We further 
constrained our sample by removing a) objects to the left (i.e. 
bluer) of the 5 Myr isochrone from Baraffe et al. dl9981 l. which 
are presumably background objects or in some cases young stars 
with strong scattered light emission from a circumstellar enve- 
lope or a nearly edge-on disk, and b) those objects located in 
the top right-hand part of the J vs I-J color-magnitude diagram, 
which are presumably foreground objects with J and I-J colors 
that cannot be explained by any, even highly extincted, 1 Myr old 
star at the distance of the ONC. Finally, our sample is reduced to 
638 periodic variables (Table 2). 

For 81 periodic variables, 64 outside and 17 inside Rciuster, we 
derived masses of < 0.075 M Q . This number is smaller than the 
124 BD candidates found in Paper 1 from I band data alone, re- 
sulting from the lower extinction in the near-infrared bands with 
the possible combination of higher variability in the I band. The 
determination of individual extinction values and spectral type 
will be necessary to unambiguously determine the masses of the 
sample. 

4. NIR excess determination 

4.1. NIR excess from (l-K) colors 

We used an index based on (I-K) colors in which I-band fluxes 
are thought to be dominated by photospheric emission, while 
K-band fluxes can have an additional strong component from 
the circumstellar disk (Hillenbrand et al. 1998). This approach 
has the advantage of a longer wavelength base, which should in 
principle result in more pronounced excess than other commonly 
used indices (e.g. A(H-K)). Unfortunately, a NIR excess derived 
from I-K colors requires the knowledge of individual extinction 
values and spectral types, both unknown for most of our ob- 
jects. As stated above, we assumed average A v values of 1.4 mag 
and 0.7 mag for objects inside and outside Rciuster, following the 
methodology stated in the previous section. We are aware that 
near the ONC center many objects may be more highly extincted 

2 As in Paper 1, we use models by Baraffe et al. (1998) for our mass 
estimations. To allow a proper comparison with previous studies, each 
time a mass value is given here it is transfered into the mass scale of the 
mentioned model. Mass values estimated from other models are given 
explicitly only if needed. 
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Table 2. Sample of the NIR catalogue of periodic variables in Paper 1. The complete catalogue is available electronically at the 
CDS. 



ID" 


Ijmag] 


*-€YY L feJ 


JTmasl 


J^ff-rmasl 


HTmael 


H^,,- rmael 

" r L ©J 


KFmasl 


K^rniael 


A(7 - 


Flag^ 


1 


15.535 


0.002 


13.540 


0.030 


12.820 


0.040 


12.510 


0.030 


-0.115 


2MASS 


3 


13.205 


0.001 


11.284 


0.001 


10.970 


0.001 


10.116 


0.000 


0.109 


UKIDSS 


7 


15.637 


0.002 


13.930 


0.030 


13.370 


0.040 


12.940 


0.020 


-0.013 


2MASS 


13 


15.551 


0.002 


14.002 


0.003 


13.374 


0.003 


13.136 


0.003 


-0.125 


UKIDSS 


10384 


17.954 


0.054 


14.960 


-1.000 


14.160 


-1.000 


13.520 


-1.000 


0.114 


VLT 


10387 


17.674 


0.011 


14.770 


0.020 


13.660 


0.030 


13.070 


0.020 


0.334 


2MASS 


20323 


14.495 


0.001 


12.260 


0.020 


11.210 


0.030 


10.810 


0.020 


0.145 


2MASS 


20325 


15.397 


0.005 


13.140 


0.040 


12.410 


0.040 


12.100 


0.040 


-0.313 


2MASS 


20326 


18.853 


0.047 


14.680 


-1.000 


13.270 


-1.000 


12.370 


-1.000 


1.683 


VLT 


20672 


19.654 


0.048 


17.276 


0.022 


15.957 


0.023 


15.638 


0.022 


0.516 


UKIDSS 



" Identification numbers from Rodriguez-Ledesma et al. (2009) 

6 I magnitudes from Rodriguez-Ledesma et al. (2009). 

c A(I - K)> 0.3 indicates the precense of a circumstellar disk. 

d This flag shows the origin of the NIR data. The possibilities are: 2MASS, UKIDSS, VLT. The VLT data set was not provided with 
individual photometric errors. A description of the estimated errors is available in Sect. 2. 



than the average, resulting in a sample of objects pretending to 
have disks, which is actually contaminated by highly extincted 
objects. On the other hand, in the outer regions a certain fraction 
of objects have lower extinction than the median value consid- 
ered here and therefore we may underestimate the fraction of ob- 
jects showing indications of a circumstellar disk in this region. 
Since we do not have information on individual spectral types, 
we need to make an assumption on the intrinsic photospheric 
color (/ - K)q. We used the observed I-J colors corrected by the 
assumed average A v values, and use the Baraffe et al. 0998) 
color information for 1 Myr old objects to derive individual in- 
trinsic (/ - K)o colors for our sample. 

We then used the definition by Hillenbrand et al. (119981 1 to ac- 
count for the amount of near-infrared excess from the I-K color: 

A(/ - K) = (I - K) measured - 0.5xA v -(I- K) , (1) 

where (/ - K)o are the intrinsic colors. 

We considered A(I-K)=0.3 mag (Hillenbrand et al. [19981 as a 
lower limit in order to be regarded as a result from a circumstel- 
lar disk. We found that about 28% of the whole sample show a 
NIR excess according to the A(I-K) index. We denoted these ob- 
jects as A(I-K)+ in Table 3. Although several assumptions were 
made to calculate the A(I-K) index, the values derived here will 
allow for the first time the comparison of the rotation-disk be- 
havior between the higher mass objects and the very low mass 
and even substellar objects. 

4.2. Extinction-free indices 

For comparison purposes we also calculated extinction-free in- 
dices. The (unknown) extinction of the individual objects and 
the IR excess from a circumstellar disk cannot be easily disen- 
tangled in color-color diagrams. As a second disk-indicator, we 
use an "extinction-free" index Q as defined by Damiani et al. 
(2007). In a J-H versus H-K diagram, there is a permitted strip 
of colors on which dwarfs should be located and which depends 
on the spectral type and the optical extinction. Since our sample 

3 To compute those indices, an extinction law has to be assumed and 
therefore they are not completely independent of extinction. In all cases 
we assumed the Rieke & Lebofsky (1985) extinction law with R= 3.1. 



includes objects in a broad mass range (i.e. from about 1.5 M to 
~ 20 Mj up ), we defined three limiting values of the index Qjhhk 
according to the three different mass regimes considered here. 
The extinction-free index provides a measurement of the devi- 
ation from these permitted strips. We then used the NIR index 
Qjhhk defined as follows 

Qjhhk = (J - H) - (H - K)(E(J - H)/E(H - K)), (2) 

in which the ratio of the color-extinctions E(J-H)/E(H-K) is 
based on the extinction law of Rieke & Lebofsky d!9851 l. 
We note that the calculated Qjhhk values provide a conservative 
limit on the presence of an IR excess and therefore the presence 
of a disk. One should keep in mind that many objects in the 
permitted strip presumably show red J-H and/or H-K colors due 
to a combination of extinction and near-infrared emission from 
a circumstellar disk (see Meyer et al. 1997). For many objects 
located at larger distances from the central Trapezium cluster, 
where the average extinction is expected to be relatively small, 
NIR colors redder than those of dwarf stars are more likely to 
be due to a circumstellar disk, and are missed by the Qjhhk- 
This conservative disk indicator explains the ~20% disk fraction 
found, a value significantly smaller than the one derived from 
the A(I-K) index. Although the Qjhhk index misses a large 
fraction of disks, we will mention the results obtained from it 
for comparison purposes. 



5. Correlation between infrared excess and 
rotational periods 

Below we investigate whether there is a correlation between 
the period distribution of the periodic variables measured in 
Paper 1 and the infrared excess. As described in the introduc- 
tion such a correlation was found in some studies for the higher 
mass periodic variables (e.g. H2002, Rebull et al. 2006, Cieza & 
Baliber, 120071) . Here we mainly want to find out whether there 
is a rotation-disk connection in the very low and substellar mass 
regime, and compare these results with higher mass objects. 
Figure[3] shows the rotational period as a function of the A(I- 
K) index for objects in the whole WFI field. Since mass effects 
may bias any existing correlation between rotation and disks we 
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Fig. 3. Rotational period as a function of A(I-K) for periodic 
variables with masses >O.4M , 0.075-0.4 M G , and <O.O75M 
for the whole WFI field. Vertical solid lines indicate A(I-K)= 0.3, 
with objects located right to this line showing an infrared excess. 
The median rotational periods for objects with and without an 
infrared excess are shown as horizontal dotted lines in all three 
mass regimes. 



perform all following analysis for the periodic variables in the 
already defined three mass bins. In the highest and intermediate 
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Fig. 4. Period distribution for objects with and without disks in 
the three mass bins investigated here (see text for details). Solid 
lines show the period distribution of objects without disks, while 
dotted lines represent the distribution for disk holders. Median 
rotational periods are shown as vertical lines. 



mass regime, Fig.|3]shows that fast rotators (i.e. P < 2 days) are 
mostly objects without NIR excess, while there is a significant 
fraction of slow rotators without indications of a disk, in addition 
to the slow rotators with NIR excesses detected. In the substel- 
lar regime no trend is clear from Fig. [3] The horizontal lines in 
each diagram indicate the median rotational periods for objects 
with and without IR excess (see Table 3). In order to give statisti- 
cal significance to the found trends we performed a two-sample 
K-S test, which shows differences in the period distribution of 
objects with and without disks in the highest mass bin at > 2 cr 
level (P = 0.02, ti\ =27, «2 = 85), in the intermediate mass bin 
at » 3 a level (P = 10~ 6 , m = 121, n z = 324), while the period 
distribution of disks and non-disks holders in the substellar mass 
bin are statistically identical (P = 0.95, n\ — 33, ni =48). 
The period distribution of the periodic variables with and with- 
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Fig. 5. The same as Fig.[3]but for the peak-to-peak (ptp) ampli- 
tude. 



out disks in the three mass bins are easier to visualise in Fig. 
[4] where the red dotted histograms represent the period distribu- 
tions of periodic variables with disks and the solid histograms 
are the corresponding period distributions for objects without 
disks. Figure|4] shows that the distribution of objects with and 
without disks in the lowest mass bins look very similar. In the 
high and intermediate mass bins objects with disks rotate on av- 
erage « 1.4-1.8 times slower than objects without disks. In the 
substellar regime no such difference in the median rotational pe- 
riods of objects with and without disks was found. 
The same trend is found when the Qjhhk index is used as diag- 
nosis of a circumstellar disk. 

In addition to mass effects, age effects may bias the resulting 
correlations, since older objects tend to rotate faster and the disk- 
fraction among an older population should also be smaller. We 
suggested a possible age spread in the ONC in Paper 1, with 
younger objects being located inside Rciuster- Therefore it would 
in principle be interesting to do the analysis shown in Fig.|4]for 
inside and outside Rciuster- However, due to the reduced number 
of objects in the two spatial regimes and, particularly in the sub- 
stellar mass regime, no statistically meaningful analysis could 
be obtained. We only show the median values of the rotational 
period distributions in Table 3, which may suggest differences 
between the two spatial regions. 
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Fig. 6. Same as Fig. [4] but for the peak-to-peak (ptp) amplitude 
distribution. Only the intermediate and substellar mass bin is 
shown since we have no information on the ptp amplitudes for 
most of the highest mass periodic variables from H2002. 



6. Correlation between infrared excess and 
peak-to-peak amplitudes 

In the following we analyze whether there is a correlation be- 
tween the peak-to-peak (ptp) amplitude of the periodic variables 
and NIR excess, i.e. the presence of a circumstellar disk. Since 
H2002 did not provide ptp amplitudes, we use here only the peri- 
odic variables measured in Paper 1 and not the combined sample 
as in previous sections (for details in how the ptp amplitudes 
were determined see Paper 1). Due to the small number of ob- 
jects with M > 0.4 M we only investigate the ptp-disk correla- 
tion in the 0.4-0.075 M Q and < 0.075 M Q mass bins. 
Figure|5] shows the ptp amplitudes as a function of the A(I-K) 
index for these two mass intervals for the whole WFI field. The 
horizontal lines in each diagram indicate the median ptp ampli- 
tudes for objects with and without IR excess. From Fig. [5] it is 
evident that objects with a NIR excess have on average about a 
factor 2-3 larger ptp amplitudes than those without a NIR excess. 
For the intermediate mass bin the two-sided K-S test gives a 10" 7 
chance that the distributions of ptp amplitudes of objects show- 
ing infrared excess and the corresponding distributions of non- 
excess objects belong to the same population («i = 93, «2 = 243). 
In the substellar mass regime the K-S test implies a difference of 
the two populations at 3cr level (P- 0.001, n\ =30, «2=43). 
Among the objects with NIR excess there seems to be a trend 
of increasing ptp amplitude with increasing amount of NIR ex- 
cess. Figure[6] helps to visualise the ptp distributions of objects 
with (dotted histograms) and without (solid histograms) disks. 
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Fig. 7. Disk fraction as a function of rotational period for the three mass regimes studied here. The disk fractions are calculated for 
rotational period bins of 2.5 days with error bars representing the statistical standard error. 



Table 3. Median values of rotational period for objects with (A(I-K) +) and without (A(I-K) -) IR excess. 





R<Rcluster 


R>Rcliistei 


Whole WFI field 




A(I-K)+ 


A(I-K)- 


A(I-K)+ 


A(I-K)- 


A(I-K)+ 


A(I-K)- 


Mass bins^ 


A B C 


A B C 


A B C 


A B C 


ABC 


A B C 


N 


13 56 8 


49 120 8 


14 65 25 


36 204 40 


27 121 33 


85 324 48 


< P >[days] 


6.3 4.7 1.5 


4.4 2.9 1.7 


6.4 3.8 1.7 


3.8 2.0 1.6 


5.8 4.5 1.7 


4.2 2.5 1.7 



" Mass bins as defined in Sect. 3: A corresponds to the bright bin containing stars with M > 0.4 M a , B to the intermediate bin with 0.4< M < 
0.075 M , and C to the faint bin in which all objects have M < 0.075 M Q . 



Table 4. Median values of ptp amplitudes for objects with (A(I-K) +) and without (A(I-K) -) IR excess. 





R<Rcluster 


R>Rcluster 


Whole WFI field 




A(I-K)+ 


A(I-K)- 


A(I-K) + 


A(I-K)- 


A(I-K)+ 


A(I-K)- 


Mass bins^ 


B C 


B C 


B C 


B C 


B C 


B C 


N 


44 8 


97 8 


48 22 


150 35 


93 30 


243 43 


< PTP >[mag] 


0.12 0.24 


0.05 0.16 


0.11 0.23 


0.03 0.08 


0.12 0.24 


0.04 0.1 



8 Mass bins as defined in Sect.3: A is not shown due to the small sample number, B corresponds to the intermediate bin with 0.4< M < 0.075 
M Q and C to the faint bin in which all objects have M < 0.075 M e . 



Median ptp values are shown in Table 4. The lack of objects 
with ptp amplitudes below 0.025 mag among the disk-holders 
is evident from both Figs. [5] and [6] We found that «79 + 2% 
of all periodic variables without disks have ptp amplitudes be- 
low 0. 1 mag, while the corresponding fraction among the ob- 
jects with disks is about 35+4%. These results reveal a very 
robust indication of differences in the variability level of objects 
with and without circumstellar disks, in which objects with disks 
show on average much larger amplitude variations than objects 
without disks. 

The same trend, median ptp values, and level of significance 
were found when the Qjhhk index is used as diagnosis of a cir- 
cumstellar disk, providing additional support for the found cor- 
relation. 

In Paper 1 we found a correlation between the rotational period 
and the level of variability for periodic variables in the ONC and 
NGC2264. As stated in Paper 1, this finding can probably be ex- 
plained by differences in the spot coverage and/or differences 
in the magnetic field topologies. The strong rotation-disk and 
ptp-disk correlations found in the intermediate mass bin studied 
here may support this scenario. We found it rather interesting 
that the correlation between ptp amplitudes and NIR excess is 



strong in the substellar mass bin, while no rotation-disk correla- 
tion was found. This finding may indicate that even dipolar-like 
magnetic field topologies, which more likely produce large spot 
groups near the magnetic poles resulting in a high level periodic 
variability, may not be strong enough to allow effective disk cou- 
pling and braking in the substellar mass regime. 

7. Disk frequency 

The disk frequency among our sample of 638 periodic variables 
can be estimated from the number of objects with a NIR excess. 
According to the A(I-K) index 28 + 2 % of them have circum- 
stellar disks. From the more conservative Qjhhk index, a lower 
disk frequency of about 20% was derived. It has to be pointed 
out that both disk indicators are biased towards a small disk 
fraction since many periodic variables are WTTSs. If all stars 
in the ONC are considered the disk-frequency is substantially 
higher, namely between 50% and 90% (Hillenbrand et al. 1998'). 
These numbers are consistent with the considerably larger frac- 
tion of irregular variables (CTTSs) relative to periodic variables 
we found in Paper 1 for the ONC (ratio of irregular to periodic 
variables « 1.7). In addition, as will be discussed in Sect.[8] the 
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derived disk fraction may also be smaller than expected since 
NIR excess is not an optimum disk indicator. 
In order to investigate a possible dependence of the disk fraction 
on the mass of the periodic variable, we derived the fraction of 
objects with and without disks in the three mass bins investigated 
here. We found that the disk frequency slightly increases from 
24% + 4% to 27% + 2% for the highest to the intermediate mass 
bins, while the increase in the disk fraction among the substellar 
objects is stronger (41% + 5%fl FigureQ shows the disk frac- 
tion as a function of rotational period, in the three mass regimes 
studied. In the high and intermediate mass bin, the disk fractions 
are computed in bins of 2.5 days up to 7.5 days, while due to 
the small number of slow rotators, the last bin includes all ob- 
jects with a rotational period larger than 7.5 days. In the substel- 
lar regime, due to the much smaller sample numbers, only three 
bins are used, i.e. P < 2.5 days, 2.5 < P < 5 days, and P > 5 days. 
From Fig. [7] is evident that there is a trend towards larger disk 
fractions for the slower rotating stellar objects, while in the sub- 
stellar mass regime no such trend is indicated at all, despite the 
fact that the disk fraction is rather high for all periods. Figure|7] 
fully supports the finding of a correlation between period and 
NIR excess for objects with masses > 0.075 M . 
In each mass regime we found that the disk frequency is always 
higher inside R c iuster than in the presumably older region outside 

Rcluster- 



8. Discussion 

8.1. Higher mass objects 

According to Sect. [3] a correlation between rotational periods 
and the presence of circumstellar disks was found in the higher 
mass regime at > 2 <x level. As outlined in the Introduction, in 
the past decade several authors have investigated correlations 
between rotation rates and the presence of circumstellar disks in 
the solar and subsolar mass regime (i.e. M > 0.2 M Q ; e.g. H2002; 
Makidon, 2004] Lamm et al.[2005 1 Rebull et al. 2006 ; Cieza and 
Baliber, 2007 ; Nguyen et al. 2009). Some studies found evidence 
for such a connection and others did not. Most of the discrepant 
results are probably related to the inhomogeneous samples that 
have been used for comparison. Due to the strong dependence 
of the rotational period with mass and age, proper comparisons 
can only be done if the samples contain periodic variables in the 
same mass regime and of similar ages. In the ONC, the mass 
effect was taken into account in the work of Cieza & Baliber 
(2007 ). They use IRAC excesses as disk indicators together with 
rotational periods from the literature, and find for objects with 
masses of ^ 0.4 M clear evidence of angular momentum regu- 
lation by circumstellar disks, confirming the correlation already 
found by H2002. 

H2002 argued for a significant correlation between A(I-K) and 
rotational periods for objects with masses of >O.2M (with 
very few objects with M < 0.2 M ) in the ONC. They tested the 
significance of the mentioned correlation by a non-parametric 
Spearman test, which resulted in less than a 10~ 10 (n = 325) 
chance that the quantities are not correlated. Since they did not 
provide the significance of the mentioned correlation for differ- 
ent mass regimes, we reanalyzed their data set and found that 
when dividing in a higher and a lower mass bin (i.e. M == 0.4 M 
or M 1 0.25 M with the D' Antona and Mazzitelli (119971 1 mod- 
els used in H2002), the level of significance of the NIR-rotation 



4 Disk fraction uncertainties correspond to the statistical standard er- 
ror [(disk fraction * (1 - diskfraction)/ N)] [/2 . 



correlation decreases substantially to 2xl0~ 3 (n= 157) for ob- 
jects with M > 0.4 M and lxlO -6 (n = 168) for the lower mass 
objects in their sample. The correlation is much stronger for the 
less massive objects in their sample, which agrees with our re- 
sults. 

The same test performed for our estimated A(I-K) index resulted 
in a 10~ 9 (n = 638), 10~ 2 (n = 112), 10~ 8 (n = 445), and 6x10" 1 
(n — 81) chance that NIR-excess from the A(I-K) index and rota- 
tional periods are not correlated for the whole sample, the higher, 
intermediate, and lower mass objects, respectively. These results 
show a trend which is consistent with H2002, in the sense that 
both higher mass objects and substellar objects tend to have a 
weaker correlation than the intermediate mass objects in our 
sample. 

We emphasise that it is not at all the purpose of this paper to 
investigate the degree of significance of a NIR excess-rotation 
correlation for higher mass stars (> 0.4 M ), but mostly to in- 
vestigate the relative change of such a correlation with mass. 
This relative change is less affected by the use of a non-ideal 
disk indicator, because the same disk indicator is used in all in- 
vestigated mass regimes. 

8.2. Substellar mass objects 

According to Figs. [3] and [4] there is no correlation between the 
NIR excess (i.e. the presence of a circumstellar disk) and the 
rotation rate in the substellar mass regime. The lack of this cor- 
relation can in principle be caused by the combination of the 
following effects: 

1 . NIR excesses are not sufficiently good accretion disk indica- 
tors. 

2. There is a certain angular momentum loss in the BD mass 
regime at the age of «1 Myr, but the following effects 
strongly weaken any correlation: 

- vast spread in rotational periods. 

- sample has a too large spread in mass and age. 

3. Angular momentum loss rates are indeed significantly lower 
in BDs than in higher mass objects. 

Below we will discuss the possible relevance of these effects in 
more detail . 

8.2.1. Are NIR excesses good accretion disks indicators? 

As outlined in the Introduction, theoretical models require on- 
going mass accretion from a circumstellar disk onto the star to 
account for angular momentum losses. To test these theoretical 
predictions observationally it is necessary to measure mass ac- 
cretion rates. Some of the more commonly used mass accretion 
indicators are the equivalent width (EW) of Ha or other suitable 
emission lines (e.g. Ca II triplet), continuum veiling, and UV ex- 
cess. Near-infrared excesses can only account for the presence of 
dust in a circumstellar disk which does not imply that gas is be- 
ing accreted from this disk onto the star. This means it is possible 
to detect a NIR excess from a purely irradiated circumstellar disk 
even if the mass accretion rate is zero (see for example Fig. 6 in 
Hillenbrand et al. 1998). Nevertheless, an impressive correlation 
between NIR excess and mass accretion indicators such as the 
Ha emission line width has been found for many YSOs over a 
broad mass range (e.g. Sicilia-Aguilar et al. 120051 12006). This 
strong correlation allows us to use NIR excess as a good accre- 
tion disk indicator in this study, since a more direct accretion 
indicator (e.g. Ha emission line flux) is not available for our 
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sample. 

Near-infrared excesses have been widely used as diagnostic of a 
circumstellar disk and have succeeded in detecting large fraction 
of very low mass stars and brown dwarfs with circumstellar disks 
(e.g. Lada et al. l2000l Muench et al. [200Tb . Muench et al. d200Tb 
found that about 50% of their brown dwarf candidates in the 
Trapezium region show NIR excesses suggestive of a circumstel- 
lar disk. This fraction is similar to the one found by Lada et al. 
(2000 ) for very low mass objects, which was increased to about 
85% when L-band (3.8 fim) observations were used. A similar 
increase in the disk fraction among the brown dwarf popula- 
tion studied here would therefore be expected from longer wave- 
length observations. In the last few years, IRAC excesses were 
used in young clusters to unambiguously determine the presence 
of a circumstellar disk (e.g. Megeath et al. 2005). The problem 
with IRAC data arises in the inner region of the ONC, where 
the nebular background and the stellar density is very high, and 
therefore a considerable number of even bright objects are miss- 
ing IRAC detections (Megeath et al. 120051 Rebull et al. [2006b . 
Up to now, only objects in the ONC with spectral types earlier 
than M2 can be precisely analyzed and searched for disk pres- 
ence and a rotational-disk connection by means of IRAC data 
(Rebull et al. 120061 and Cieza & Baliber l2007l l. We then have to 
rely on the less favourable NIR excesses as indicators of disks 
around very low mass and substellar objects. 
We like to remind the reader that the photometric method used 
for measuring periodic brightness modulations is more efficient 
among the WTTSs than for CTTSs, since irregular variations 
can add a strong noise to the periodic modulation, preventing 
the detection of the periodic signal. Due to this noise we are def- 
initely missing periodic variables among the young highly ac- 
tive CTTSs, which are expected to have larger amounts of NIR 
excess (i.e. larger mass accretion rates) than the more common 
WTTSs in our sample of periodic variables. This bias is an ad- 
ditional factor which prevents us from finding a stronger corre- 
lation between rotational periods and the presence of a circum- 
stellar disk. 

In summary, we believe that the used NIR disk indicator is not 
a major cause for the lack of a disk-rotation correlation found in 
the substellar regime. 

8.2.2. Correlation weakening/destroying effects 

Below we discuss the effects that can weaken or even destroy a 
correlation between rotation and mass accretion, even in the case 
of an "ideal" mass accretion indicator. 



(i.e. the <x values are a factor « 1.6 larger). In addition, even 
when the distribution is broader, the number of outlier^ 
is higher in the BD mass regime than in the intermediate 
mass regime, with 14% and 3% of outliers in each region, 
respectively. 

The larger spread in rotational periods among the BDs in 
our sample seems to be an important correlation decreasing 
factor in the BD regime. 

- Mass and age effects: 

As stated in Sect. 8.1, the correlation between rotational pe- 
riod and NIR excess can be biased by the strong dependence 
of rotation with mass and age. 

If a sample of periodic variables with a broad mass range 
is used in order to find a correlation with NIR excess the 
following effect may mimic a rotation-disk connection: 
higher mass objects rotate on average slower and show on 
average larger NIR excesses than lower mass objects, the 
latter due to the higher contrast between the photosphere 
and the inner disk emission (e.g. Littlefair et al. 2005 ). 
An age spread may also mimic an intrinsic correlation 
between rotational periods and disks. 

The intermediate mass regime has a mass range of about 
0.4/0.075 ~ 5. However, the mass range in the BD regime 
is most likely much smaller since we have only very few 
objects below 0.025 M , i.e. the mass range is there about 
3 (0.075/0.025). Due to the smaller mass range in the BD 
regime it cannot be excluded that this effect contributes to 
the lack of a correlation between disks and rotation rates. To 
test this hypothesis, we divided the intermediate mass regime 
in two mass bins, namely objects with 0.4<M< O.17M 
(intermediate-high) and 0.17<M < 0.075 M Q (intermediate- 
low), and analyzed the new period distributions of objects 
with and without NIR excess. The KS test resulted in a 99% 
and 92% probability that the period distributions of objects 
with and without NIR excess belong to different populations 
in the intermediate-high and intermediate-low mass bins, 
respectively. The mass ranges in both these newly defined 
mass bins are about 2.3 and therefore smaller than the mass 
range in the substellar regime; nevertheless, rotation-disk 
correlations are found to be more significant than in the 
substellar mass regime. The statistical significance of the 
mentioned correlation decreases toward lower masses, 
which most probably indicates that mass effects are not 
relevant within the mass intervals studied here. 



- Large spread of rotational periods: 

From Fig. [3] it is clear that the period distribution of the 
periodic variables in all three mass bins show a vast spread. 
If this spread can weaken or even destroy a correlation 
between rotation and NIR excess, it should in principle 
affect all three mass bins. To test whether there is a larger 
spread in rotational periods in the BD regime than in the 
intermediate mass regime, in which a highly significant 
correlation was found, we computed P/Pmedian values and 
analyzed the median and standard deviation (cr) of these 
P/Pmedian distributions. In addition, a log-normal distribu- 
tion was also fitted and both the geometric medians and the 
multiplicative cr were compared. The cr values provide the 
width of the distribution and therefore give information on 
the spread of the values. We found that the distribution of 
P/Pmedian in the substellar mass regime is broader than the 
corresponding distribution of the intermediate mass objects 



8.2.3. Are angular momentum loss rates smaller in BDs than 
in higher mass objects? 

We discuss here possible physical processes which can account 
for smaller angular momentum loss rates in the substellar mass 
regime compared to higher mass objects. 

- Low ionization: 

All hydro-magnetically driven angular momentum loss 
mechanisms require a sufficiently high ionization of the disk 
and/or wind involved in the angular momentum transport 
from the star and/or disk to this circumstellar matter. 
Without a sufficient ionization the magnetic field would only 

5 Number of objects with values outside the following range: Q r 
1.5x(Q 3 -Q!) and Q 3 + 1.5x(Q 3 -Q 1 ), where Qi and Q 3 are the first and 
third quartiles. For our sample, all outliers have values higher than the 
upper limit, i.e. longer periods. 
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poorly couple to the circumstellar matter, and therefore mass 
accretion rates and angular momentum transport would be 
correspondingly smaller. Since the X ray and UV fluxes 
of BDs relative to their bolometric fluxes (i.e. L x /Lb i) are 
certainly much lower than for higher mass stars, it is quite 
conceivable that these effects reduce the angular momentum 
braking in the BD regime. This low ionization/magnetic 
coupling consideration can presumably explain the strong 
decrease of M with decreasing M star discussed in the next 
section. 

- M versus M star relation : 

Observations over a broad mass range from intermediate 
mass stars down to brown dwarfs reveal a very strong 
dependence of the mass accretion rate on stellar mass with 
the global trend M<xM 2 slar (Hartmann et al. 2006]). An 
even stronger dependence, MocM 2 ^ has been reported in 
the subsolar mass regime down to the substellar limit by 
Fang et al. (2009). This dependence explains the very small 
mass accretion rates of < 1O~ 1O M measured in very low 
mass and substellar objects (e.g. Muzerolle et al. 2005). 
Obviously this non linear relation between mass accretion 
and stellar mass indicates that the relative braking J/ J 
is strongly mass dependent, with a substantially smaller 
braking in the lower mass objects. 

In addition, Hartmann et al. (2006) have shown that dM/dt 
is probably also highly mass dependent in the sense that it 
strongly increases with decreasing mass; i.e. M decreases 
much faster with evolutionary age in BDs than in solar mass 
stars. 

- Relevant time scales for disk-locking: Ti oc k and To'- 

If one assumes that the disk locking scenario explains the 
rotational data one can consider two extreme cases: case 1; 
all objects are locked to their disks at their birth, and the 
amount of angular momentum loss strongly depends on how 
long this locking is effective (i.e. how long T/ oc t) is, and case 
2; they are not locked to their disks at early stages and it takes 
a certain time To until they are locked. 
In the first case, the above discussed M/dt versus mass con- 
sideration do indeed argue that t/ oc a in BDs is much shorter 
than in solar-mass stars, since in the former case M will de- 
crease much more rapidly. Therefore BDs are expected to 
be braked for a shorter time than higher mass objects, after 
which they spin up conserving angular momentum or losing 
angular momentum in a moderate way. 
Concerning case 2, Hartmann ([20021 argued that disk- 
locking is not instantaneous but a certain time To is needed 
before disk-locking is achieved. Equation [3] (which is Eq. 8 
in Hartmann 2002 ) shows that this time is proportional to the 
mass, the angular velocity, and to the inverse of the accretion 
rate: 

t d £ 4.5xW 6 yrf-^-, (3) 

in which M0.5 and M~ 8 are the mass and mass accretion 
rates in units of 0.5 M and 10~ 8 M yr~', respectively and/ 
is the angular velocity in units of the breakup velocity. 
Typical values of these quantities for a 0.5 M star are 
M^lO^Moyr 1 and f =0.1 (see Herbst et al. l200TT) . which 
result in a td* 0.5 Myr. Since M in BDs is highly time 
dependent To will strongly depend on the evolutionary 
phase. If one assumes an optically observable » 1 Myr 
old BD, a typical M = 5xlO" lo M yr _1 is expected. In this 



case To is about 10 Myr (for f«0.7 which is larger in BDs 
because they rotate close to break up velocity). This all 
means that disk coupling in BDs can only be achieved 
within in a reasonable time scale in very early evolutionary 
phases, when M was presumably two to three orders of 
magnitude higher. 



8.2.4. Concluding remarks 

There is strong observational evidence that the mass accretion 
rate M for low mass stars and substellar objects is highly mass 
dependent (M oc M 2 W 2 - 8 ), and therefore one would expect that 
the ratio // J depends on mass in a similar manner. With these 
dependencies it is not surprising that no correlation between the 
NIR excess and the rotation rate was found in the substellar 
regime. 

Before reliable conclusions can be drawn, independent disk in- 
dicators (e.g. IR excess at longer wavelengths) should be cor- 
related with rotation and, if possible, future substellar samples 
should include objects at slightly different ages. 

9. Summary 

We correlated rotational periods and peak-to-peak (ptp) ampli- 
tudes with NIR excess in a large sample of VLM stars and BDs 
in the ONC. Our results are summarized in the following: 

1 . We analyzed 638 periodic variables in the ONC from Paper 1 
and H2002 with IJHK measurements. 

2. We estimated masses from the J versus I-J color-magnitude 
diagram and found 81 objects with masses in the substellar 
mass regime. 

3 . We calculated the amount of NIR excess from a circumstellar 
disk by means of the A(I-K) index and and, for comparison 
purposes, we also computed an extinction-free index Qjhhk- 

4. We performed all our analysis in three mass bins, namely ob- 
jects with masses >O.4M , 0.4-0.075 M , and < 0.075 M . 
Our results show the existence of a highly significant 
rotation-disk correlation for the intermediate mass objects, 
in which objects with NIR excess tend to rotate slower than 
objects without NIR excess. In the highest mass bin there is 
a correlation at > 2 cr level. Interestingly, and most impor- 
tant for this study, we found no correlation in the substellar 
regime. 

5. We found a highly significant correlation between the ptp 
amplitudes and the NIR excess in the two mass bins (0.4- 
0.075 M @ , and < 0.075 M ) with available ptp values, in 
which objects with circumstellar disks show on average a 
factor two to three larger ptp amplitudes than those without 
indications of a disk. 

6. We estimated a disk frequency between 20% and 28% 
among the 638 periodic variables in the ONC. This fraction 
is smaller than the one estimated for all objects in the ONC 
by Hillenbrand et al. d 19981 1. which is expected because (a) 
the photometric method used to derived rotational periods is 
more effective among the WTTSs than for the CTTSs, with 
the former group having a much smaller disk frequency, (b) 
NIR excesses miss a certain fraction of objects with circum- 
stellar disks. 

7. The three mass regimes studied show different rotation-disk 
behaviors which may be attributed to the different physical 
properties of objects in each mass bin. 

The first division line (0.4 M ) coincides with the onset of 
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full convection. It may be possible that the different be- 
havior between high and intermediate mass objects in our 
sample are related to the different ways magnetic fields are 
produced, resulting in different topologies, spots coverages 
and/or magnetic field strengths. 

The second boundary is the substellar limit. The strong de- 
pendence of mass accretion rate on mass seems to be crucial 
to understand why the lower mass objects in our sample suf- 
fer smaller angular momentum losses than those with larger 
masses. Our observations suggest that at an age of about 
1 Myr disk-braking does not influence the rotational rate of 
most substellar objects. 

In order to exclude possible observational biases by the 
method used here, future studies should consider IR excess 
at longer wavelengths or a direct mass accretion indicator. 

Acknowledgements. The authors thank Mark MacCaughrean for providing the 
VLT-ISAAC NIR data used for 1 13 objects. We especially thank the anonymous 
referee for her/his valuable comments on the original version of this manuscript. 
M.V.R.L. acknowledges support of the International Max Planck Reseach School 
for Astronomy and Cosmic Physics of the University of Heidelberg. 
This publication makes use of data products from the Two Micron All Sky 
Survey, which is a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/California Institute of Technology, 
funded by the National Aeronautics and Space Administration and the National 
Science Foundation. This work is based in part on data obtained as part of the 
UKIRT Infrared Deep Sky Survey. 



Makidon, R. B., Rebull, L. M., Strom, S. E., Adams, M. T., Patten, B. M. 2004, 

AJ, 127, 2228. 
Matt, S., & Pudritz, R., E. 2005 ApJ 632 L135-L138. 

Megeath, S.T., Flaherty, K.M., Hora, J., Allen, L.E., Fazio G.G., et al. 2005, 
IAUS 227, 383. 

Meyer, M. R., Calvet, N., & Hillenbrand, L. A. 1997, AJ, 114, 288 
Mohanty, S., Jayawardhana, R., Basri, G. 2004 ApJ, 626, 498. 
Muench, A. A., Alves, J., Lada, C. J., Lada, E. A. 2001, ApJ, 558L, 51. 
Nguyen D. C, Jayawardhana R., van Kerkwijk M., Brandeker A., Scholz A., et 

al. 2009 ApJ 695,1648. 
Rebull. L. M. 2001, AJ, 121, 1676. 

Rebull. L. M., Stauffer, J. R., Megeath, T., Hora, J., Hartmann, L. 2006, ApJ 
646,297. 

Rieke, G. H., Lebofsky, M. J 1985 ApJ, 288, 618. 

Rodrfguez-Ledesma, M. V., Mundt R., Eisloffel, J. 2009 A&A, 502, 883. 
Scholz, A., & Eisloffel, J. 2004 A&A, 419, 249. 
Scholz, A., & Eisloffel, J. 2005 A&A, 429, 1007. 

Sicilia-Aguilar, A., Hartmann, L. W., Szentgyorgyi, A. H.,et al. 2005 AJ, 129, 
363. 

Sicilia-Aguilar, A., Hartmann, L. W., Frsz, G., Henning, T, Dullemond, C, et 

al. 2006 AJ, 132,2 135 
Shu, F, Najita, J., Ostriker, E., Wilkin, E, Ruden, S. et al., 1994 ApJ, 429, 781. 
Vogel, S. N., Kuhi, L. V. 1981, ApJ, 245,960. 



References 

Baraffe, I., Chabrier, G., Allard, E, Hauschildt, P. H. 1998 A&A, 337, 403. 
Bessell, M. S. & Brett, J. M. 1988 PASP, 100, 1134. 

Biazzo, K., Melo, C. H. E, Pasquini, L., Randich, S., Bouvier, J. et al. 2009 

2009arXiv0909.5035B. 
Bodenheimer, P. Ann. Rev. A&A, 33, 199-238 
Bouvier, J. 2007 IAUS, 243, 23 IB 
Camenzind, M., 1990 RvMA, 3, 234. 
Carpenter, J. 2001, AJ, 121, 2851. 

Casali, M., Adamson, A., Alves de Oliveira, C, Almaini, O., Burch, K., et al. 

2007, A&A 467, 777C. 
Chaisson, E. J., Dopita, M. A. 1977, A&A, 56, 385. 
Cieza, L., Baliber, N. 2007, ApJ, 671, 605 

Damiani, E, Prisinzano, L., Micela, G., Sciortino, S. 2006, A&A, 459, 477. 
D'Antona, E, Mazzitelli, I., 1997, in Cool stars in clusters and associations, ed. 

R. Pallavicini, & G. Micela, Mem. Soc. Astron. Ital.,68 
Edwards, S., Strom, S. E., Hartigan, P., et al. 1993 AJ, 106, 372. 
Ghosh, P., & Lamb, F. K. 1978, ApJ, 223, L83 

Fang, M., van Boekel, R., Wang, W., Carmona, A., Sicilia-Aguilar,A., et al. 2009 
A&A 504, 461 

Ferreira, J., Dougados, C, Cabrit, S. 2006 A&A, 453, 785F. 
Hartmann, L. 2002 ApJ, 578, 914. 

Hartmann, L., D'Alessio, P., Calvet, N., Muzerolle, J. 2006 ApJ, 648, 484H. 

Herbst, W., Bailer-Jones, C.A.L., Mundt R. 2001 ApJ, 554, L197. 

Herbst, W., Bailer- Jones, C.A.L. Mundt R., Meisenheimer, K., Wackermann, R. 

2002 ApJ, 396,513. 
Herbst, W., Eisloffel, J., Mundt R. Scholz A., 2007 prpl.conf, 297. 
Hewett, P. C, Warren, S. J., Leggett, S. K., Hodgkin, S. T. 2006, MNRAS, 367, 

454 

Hillenbrand, L. 1997 AJ, 113, 1733. 

Hillenbrand, L., Strom, S., Calvet, N., Merrill, K., Gatley, I., et al 1998, AJ, 
116,1816. 

Kirkpatrick, J., Reid, I., Liebert, J., Gizis, J., Burgasser, A. et al. 2000, AJ, 120, 
447. 

Konigl, A. 1991 ApJ, 370, L39. 

Lada, C. J., Muench, A. A., Haisch, K. E. Jr., Lada, E. A., Alves, J. E, et al. 2000 
AJ, 120.3162. 

Lamm M., Bailer-Jones, C.A.L., Mundt R., Herbst, W., Scholz A. 2004 A&A, 
417,557. 

Lamm M., Bailer-Jones, C.A.L., Mundt R., Herbst, W. 2005 A&A, 430, 1005. 
Lawrence, A., Warren, S. J., Almaini, O., Edge, A. C, Hambly, N. C, et al. 2007, 

MNRAS,379,1599L. 
Littlefair, S. P., Naylor, T, Burningham, B., Jeffries, R. D. 2005, MNRAS, 

358,341. 



1 

0.75 

0.5 

0.25 


1 

0.75 
0.5 



M>0.4 MQ 




1 

0.75 
0.5 
0.25 




0.075<]Vy0.4 M© 



0.25 m m 



i i r 

M<O.O75JVI0 



J I L 



4 8 12 
P[days] 



